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Abstract 

We study a family of polynomials that are orthogonal with respect to the 
weight function e‘““^ in [—1,1], where cj > 0. Since this weight function is 
complex-valued and, for large uj, highly oscillatory, many results in the classical 
theory of orthogonal polynomials do not apply. In particular, the polynomi¬ 
als need not exist for all values of the parameter lo, and, once they do, their 
roots lie in the complex plane. Our results are based on analysing the Hankel 
determinants of these polynomials, reformulated in terms of high-dimensional 
oscillatory integrals which are amenable to asymptotic analysis. This analysis 
yields existence of the even-degree polynomials for large values of cj, an asymp¬ 
totic expansion of the polynomials in terms of rescaled Laguerre polynomials 
near ±1 and a description of the intricate structure of the roots of the Hankel 
determinants in the complex plane. 

This work is motivated by the design of efficient quadrature schemes for 
highly oscillatory integrals. 


1 Introduction 

Polynomials orthogonal with respect to the weight function e“^^ in [—1,1] have been 
proposed in (Asheim, Deafio, Huybrechs & Wang 2014) as a means to derive complex 
Gaussian quadrature rules for highly oscillatory integrals. Though a wealth of obser¬ 
vations were made regarding these polynomials, the focus in (Asheim et al. 2014) is on 
the numerical analysis of the quadrature error for such integrals. The current paper 
revolves around the analysis of the polynomials themselves, in the process confirming 
rigorously several of the observations that were made before. 

We proceed by analysing the Hankel determinants associated with the polynomi¬ 
als as a function of the parameter lo . Though Hankel determinants are usually men¬ 
tioned in the classical references for the theory of orthogonal polynomials, for instance 
(Szegd 1939), (Gautschi 2004), (Ghihara 1978), (Ismail 2005), and they are of inter¬ 
est in connection with integrable systems and random matrix theory, they are rarely 
a starting point for further analysis: insofar as classical orthogonal polynomials are 
concerned, it follows from standard results on the Hamburger moment problem that 
Hankel determinants are always positive and they throw little added insight on the 
underlying problem. The main reason for pursuing Hankel determinants in this paper 
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is an old result of Heine (Ismail 2005, §2.1) that they can be written as a multivariate 
integral, as shown in which in the context of our oscillatory weight function be¬ 
comes an oscillatory integral. The integral can be expanded asymptotically in negative 
powers of w and this expansion is very revealing. Furthermore, several quantities of 
interest related to the orthogonal polynomials, including the polynomials themselves, 
can be formulated in terms of Hankel-like determinants. 

To establish our notation, we consider monic orthogonal polynomials (OPs) Pni^) 
or Pn{-), formally defined as follows: 

J p‘:^{x)x’^e^‘^^ dx = 0, fc = 0,1,..., n — 1, (1-1) 


where n € and a; > 0. For simplicity of notation, in the sequel we omit the 
parameter uj and write (cc) directly. 

The standard construction of orthogonal polynomials with respect to a measure p 
involves the Hankel matrices built from the moments of the measure. In the present 
case, dp{x) = e“^^ dx, let 


/r„(w) = / dx, n G 


be the moments of the weight function and set 


Hr,. = 


ho 

hi 

l^n 

hi 

h2 

Mn+l 


hn+1 

M2n 


and 


hn = detHn, n G N, (1.2) 


the nth Hankel matrix and determinant, respectively. Note that both and hn 
depend on w, but we omit this for simplicity of notation. It is well known that the 
polynomial p„ {x) itself can also be written in terms of determinants as 



ho 

hi 

' f^n—1 

1 

—— det 

hi 

h2 

' l^n 

X 

hn-1 

l^n 

Mn+1 

h2n-l 



see for instance (Ismail 2005, Chapter 2). Note that, given a value of w, Pn{x) exists if 
and only if hn-i{uj) ^ 0. This is guaranteed whenever /r is a positive Borel measure, 
so for classical polynomials is always strictly positive, and in that case any deeper 
study of the Hankel determinants is usually of little further consequence. In our 
setting, however, the study of the Hankel determinants, and in particular the study 
of their roots, is extremely insightful. 

Note that, despite the fact that the weight function is complex-valued, the deter¬ 
minants hn are real for all n G Z_|_. This is so because 

p„{-uj) = (-l)"^„(w), n G Z+. 
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Thus, if we consider the Hankel matrix Hn{—uj), we can pull out a factor (—1)^ from 
the fc-th row and (—1)^ from the £-th column. As a consequence, 

hn{uj) = hn{-uj) = = hn{oj), 

and the result follows. 

In §2 we sketch some properties of our orthogonal polynomials and in §3 we express 
their Hankel determinant as a multivariate integral. The heart of our investigation is 
§4, where we analyse Hankel determinants and their asymptotic expansion for w ^ 1. 
This is used in §5 to explain the kissing pattern of zeros, in §6 to investigate the highly 
nontrivial behaviour of orthogonal polynomials near the endpoints ±1 and in §7 to 
shed light on the complex zeros of Hankel determinants. Finally, in §8 we extend the 
important fact that the existence and uniqueness of is assured for w > 0 small 
enough or large enough to all w > 0. 

The origins of ‘kissing polynomials’ are in highly oscillatory quadrature but, we 
believe, their behaviour and the crucial differences between them and conventional 
orthogonal polynomials are a matter of an independent mathematical interest. In this 
paper we consider the simplest possible case but we expect to return to this theme in 
future papers, investigating polynomials orthogonal with respect to the more general 
sesquilinear form, 

{fij2) = j /i(x)/2(a;)w(a;)e’‘^®("’) dx, 

where w and g are given functions of suitable regularity. Preliminary results indicate 
a wealth of fascinating and beautiful behaviour. 


2 Properties of the orthogonal polynomials 

We recall few interesting properties and observations on the polynomials p„ from 
(Asheim et al. 2014). Firstly, they are symmetric with respect to the imaginary axis: 

Pn{z) = (-l)>„(-z), z e C. (2.4) 

Note that the map z —>■ —z represents a reflection with respect to the pure imaginary 
axis: if z = x + ip then —z = —x + ip. 

Additionally, since the orthogonality relation (HtD is non-Hermitian, the monic 
OPs satisfy a three term recurrence relation, 

Pn+l(x) = (x - an)pn{,x) - /3„p„_i(x), (2.5) 

provided that three consecutive polynomials exist for a given value of lo. The initial 
values are taken as p_i(x) = 0, po{x) = 1 and the coefficients a„ and /3„ are complex 
valued. 

The coefficients of the recurrence relation can be given in terms of Hankel deter¬ 
minants, namely. 





hn — 1 J 


Pn = 


hnhn—2 

K .-1 


( 2 . 6 ) 
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where h'^ indicates differentiation with respect to uj. 

As a function of the parameter w, the recurrence coefficients themselves satisfy the 
differential-difference equations 

a; = i(/3„+i-/3„) (2.7) 

Pn — i/^n(Q^n 

This is nothing but a complex case of the classical Toda lattice equation subject 
to the Flaschka transformation, which is known to govern the deformation of the 
recurrence coefficients whenever the measure of orthogonality is a perturbation of a 
classical one with an exponential factor linear in the parameter, which is w in our case. 
We refer the reader to (Ismail 2005, §2.8) for more details. 

Another important consequence of the weight function not being positive is the fact 
that, even when Pn{x) exists for some values of n and w, its roots lie in the complex 
plane. They come in pairs of two, symmetric with respect to the imaginary axis, as a 
consequence of (12.41) . 

When u) = 0, pn is a multiple of the classical Legendre polynomial and the roots are 
real and inside the interval [—1,1]. For increasing values of ui, they follow a trajectory 
in the upper half plane illustrated in Figure HT] The trajectories corresponding to 
polynomials of consecutive even and odd degree touch at a discrete set of frequencies 
uj: the polynomials ‘kiss’, see Figure [221 



Figure 2.1: The zeros of pe (dark) and of P 7 (light). 

Regarding the asymptotic behaviour of the polynomials p„, it is of interest to con¬ 
sider the case when one of both parameters, n and w, are large. The large n asymp¬ 
totics can be deduced using for instance the Deift-Zhou steepest descent method 
applied to the corresponding Riemann-Hilbert problem, as in (Deano 2014) and ref¬ 
erences therein, see also the general monograph (Deift 2000). This analysis focuses in 
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Figure 2.2: Close-ups of the ‘kissing’ patterns near the right endpoint -1-1. 


the case when n is large and lo is fixed or grows linearly with n, and provides exis¬ 
tence oipn{z) for large enough n together with asymptotic behaviour for z in different 
regions of the complex plane. 

The asymptotic analysis for fixed n and large w is much less standard. It is 
conjectured in (Asheim et al. 2014) that the even-degree polynomials asymptotically 
behave like a product of Laguerre polynomials centered around the endpoints ±1: 


P2n(a;)-"f—j L„(-iw(x -I- l))L„(-ia;(a:: - 1)), ui ^ oo, (2.8) 


where L„ is the nth Laguerre polynomial with parameter a = 0, see (Szegd 1939). 
Thus, for large w, it seems that the orthogonal polynomials of even degree become 
approximately a product of lower degree orthogonal polynomials. This conjecture also 
implies that the roots shown in Figure ?? behave like ±1 — where c is a root of 
the Laguerre polynomial L„- 
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3 Properties of Hankel determinants 


3.1 An expression as a multivariate integral 

We commence by revisiting an old result of Heine (Ismail 2005, Section 2.1), (Szego 
1939, pg. 27) which is fundamental to our analysis. Equally fundamental is the 
method of proof, which will be reused and generalised in the sequel. 


Lemma 1 For every n £ Z+ it is true that 

0<fc<.£<n—1 


^n—1 — 


\ t t t n da;„_i • ■ • dxi dxQ. 

■ J-lJ-l -^-1 o<fe<£<n-l 

(3.9) 

Proof We write the determinant in the following form: 


hn-i = det 


/-I 

r*l 


dxo 
dxo 


dxi 


'-1 

/•i 


'-1 

/*! 


dxi 


n-l AujXn-i 
^n-1^ 


n l(jJXn — l 


Xn-ie 


da;„_i 


da;„_i 


/-I 


f-1 


f-1 




dxi 


det 




-ij-i j-i 


n—1 

n 


fc =0 


Xi 

^n-1 


a-i 

•^n-1 


a-i 

2n-2 

‘^n-1 J 

giu;® 1 

1 

1 

1 

Xo 

Xi 

^n—1 

n — 1 

^n—1 

^n—1 

Xq 

Xi 

^n-l 


^2n-2^iujXrt-i , 


■dx„_i• • • dxi dxo 


^ dx„_i• • • dxi dxo 


j-I 


\ n—1 

n 

fe=0 0<k<i<n-l 


W_^k It (a;^-a:fc)e“^“=^^dx„_i dxidxo. 


using the well known formula for the determinant of a Vandermonde matrix, and 
the notation Ax is used for the linear function xq + xi + ... + x„_i. Let tt be a 
permutation of (0,1,..., n — 1). Then, changing the order of integration, 


/ I p 1 n-l 

/ It 

-1 -^-1 k=0 


A 

^Tv{k) 


It (a^£ aj/„)e dx^_i • • • dxi dxo, 

0<k<e<n-l 


where (j{7z) is the sign of the permutation. Averaging over all n! permutations, 

hn—\ — i" f '' I ; a^n—i) 11 (a^£ ^k')^ dxyi—i • • • dxi dxi 

n. y_i o<k<e<n-i 


-1 
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where 


n —1 

g{xQ,...,Xn-i)= XI n 

TTGlIn k—0 

and n„ is the set of all the permutations of (0,1,..., n — 1). And now comes the 
serendipitous step that we will use time and again in the sequel: the observation that 
g is a determinant - specifically, a determinant of an n x n Vandermonde matrix. 
Therefore 

g(xo,...,Xn-i) = (xt-Xk) (3.10) 

0<k<e<n-l 

and the proof of ()3.9I1 is complete. □ 

3.2 The asymptotic expansion of hn 

The oscillatory integral 13.91 can be expanded asymptotically in inverse powers of w. 
Indeed, the integrand has the canonical form of a non-oscillatory function 

F{x) multiplying an oscillatory exponential, with the so-called oscillator g{x) - in this 
case simply the linear function a;o-l-a:i-|-.. ■ + Xn-i- It is well known how to derive such 
expansion, for example using repeated integration by parts (see, e.g., (Wong 2001)). 
This is straightforward in principle, but hampered by lengthy algebraic manipulations 
in our high-dimensional setting since [3)9] is an n-fold integral. In the following, we will 
use the multi-index notation of (Iserles & Nprsett 2006) to control the complexity: 

[ Fix)F‘^^^^dx (3.11) 

-'[- 1 , 1 ]" 

OO ^ 

~ (-1)” V TZiuwsr V V 

m=0 ' \k\=mvGVn 

We let V„ be the set of the 2” vertices of the n-cube [—1,1]". The index function s{v) 
of 1 ) € Vn is the number of —1 therein. For example, for n = 3 the index of the vertex 
(—1, —1,1) would be 2. 

Note that each term in the expansion, corresponding to some negative power of 
w, consists of summing over all partial derivatives of a certain total order m over all 
possible vertices of the cube [—1,1]". One may think of these derivatives as originating 
from the integration by parts technique, and they are evaluated at the vertices because 
the endpoints of all univariate integrals involved are either -|-1 or —1 and in this case 
the integrand has no singularities or stationary points. 

Specializing (13.111) to the integral (13.9p at hand, and using (I3.10|) . we note that the 
non-oscillatory function 


A TT {xt-Xkf, neZ+. (3.12) 

0<fe<^<ra-l 

is a polynomial of total degree (n—l)n. This implies that expansion (13.1111 terminates, 
as all derivatives vanish once m > in — l)n -|- 1. Since the expansion of hn-i starts 
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with oj " 


because of Lemma [TJ we expect it to have the form 


hn—l 


E 




(3.13) 


The reason for the upper bound is that for the last significant value of = (n — l)n, 
we have i + n = n'^. 

By direct calculation we find the first few expansions 


hoiio) = 


2 sinw 


w 


, , , 4 2(cos2w—1) 

hii^) = - + - -4- 

32 since 64cosa; 96 since 32sin^ce 

/i2(ce) = -g-g- h-=- 5 - , 




256 512(cos2ce-4) 3072sin2ce 768(11 cos2ce - 2) 9216sin2ce 


,10 


,11 


,12 


+ 


,13 


6912(cos2ce “ 1) _j_ 576(cos2ce — 1)^ 


ce 


14 


ce 


16 


The upper bound in (13.131) is sharp. However, the leading powers in these expressions 
are substantially higher than predicted by (13.131) and the discrepancy becomes more 
pronounced as n increases. Instead of (13.131) we have 


hn — l — ^ ^ 


hn — 1 ,. 


e=s„ 


UJ’- 


(3.14) 


where 

= 1, 52 = 2, ^3 = 5, ^4 = 8, <55 = 13, ^6 = 18, <57 = 25. 


Note that for even n the values h^-i^Sn in tbe examples above are positive constants: 
this implies the existence of even-degree polynomials for sufficiently large uj. The 
factor sincj, however, appearing for odd n, demonstrates asymptotic non-existence of 
odd-degree polynomials, approximately at integer multiples of tt. 

In view of the fact that expansion ()3.14|) terminates, it is actually an exact and 
explicit expression for hn-i- Though it may seem to have a high-order pole at oj = 0, 
the singularity is removable, since from (13.91) it follows that /i„_i is an analytic function 
of LO. 

The case Ct; = 0 corresponds to the Legendre weight function and the corresponding 
Hankel determinants are all positive. Actually, as oj —>■ O'*", we have 


ho{uj) 

hi{uj) 

h2{ui) 

h^iu)) 




- OJ 

3 45 

32 16 


315 


-UJ 


0{uj^ 


+ -—cj^ + 0(cc®), 


256 


23625 1488375 


2025 
2048 2 


1024 


UJ 


11694375 


UJ 


0{u^). 




























The quest for the leading order term in expansion (13.141) revolves around the study 
of the derivatives of the integrand at the vertices of the hypercube [—1,1]". In partic¬ 
ular, it is clear from the explicit expression (13.121) that the integrand vanishes to some 
order whenever two coordinates xt and Xk coincide. This is the case at the vertices 
and, loosely speaking, the order is determined mostly by the difference between the 
number of -I- Is and —Is at the vertex. 


4 Asymptotic analysis of a symmetric integral 

It is convenient for later developments to aim for a slightly higher level of generality 
than in 43.21 before simplifying again to the case at hand. We set 

F{x) = f{x)g{xf, (4.15) 

where the function g is the Vandermonde determinant: 

n—1 

g{x) = g{xo,xi,... ,Xn-i) = H 11 

0<k<e<n-l TT k=0 

and 7r € n„ is the set of all permutations of length n, acting for example on the 
n-tuple (xo,..., Xn-i)- In the sequel we will also assume that the smooth function / 
is symmetric in its arguments. 

Thus, we set out to study the n-fold integral 

In[f] = —^f ■ ■ ■ f f (xo,... dxo • ■ • dx„_i, (4.17) 

f (^0; ■ • ■ ; —1) ^ dxo * ' * cIXti, — ! , 

0<fc<^<n—1 

The cases of primary interest in this paper are the following: 

• / = 1, since /n[l] corresponds precisely to according to the integral (13.9|) . 

• fix) = nm=o( X — Xm), that corresponds to the polynomial hn-ip^ix). This 
follows from Heine’s formula for the orthogonal polynomial: 

1 pi pi n-l 

K(^) = ^t - / ■••/ n dxo 

^'■hn-lJ-1 0<k<An-l 

see, e.g., (Ismail 2005, Theorem 2.1.2). 

Of course, once / is a polynomial, so is F, and in that case (13.111) terminates 
and is no longer an asymptotic expansion but an exact formula. In any case, we are 
faced with computing derivatives of F(x) and, in particular, the first non-vanishing 
derivative of f at a vertex of the hypercube. 


dx^^—i, 


(4.18) 
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4.1 Vanishing derivatives of the function F 

Everything is symmetric, hence we can assume that for a vertex that contains r com¬ 
ponents equal to —1, i.e. s{v) = r, we have 

r times n—r times 

« = (4.19) 

Let us consider the following factorization: 

r—1n—1 

ar{x)= {xi-Xk), and PnA^) = 'WW_i.xi - Xk), (4.20) 

0<A:</<r—1 A:—0 l—r 

then 

F(^x) = f {Fj ex (Xq^ ■ • ■ 1 X'r—l)cen—r(Xr: • ■ ■ : ^n —l)/5n,r (^0 ; . . . ; Xn-l)- (4.21) 

We know by construction that does not vanish at a vertex, in fact j3n,r{v) = 
2 r{n-r) ^ We assume, at least for the time being, that f{v) does not vanish either, a 
condition that is clearly satisfied in the cases enumerated before. Therefore, we need 
be concerned just with at —1 and Un-r at -1-1. Since a{x + cl) = a(x) for all 
X G K", c € K, because a only depends on the differences between elements of a;, it is 
sufficient to examine these expansions at a; = 0. 

By the definition of a determinant, 

ar(x) = VDM(a:o,.. ■,Xr-i) = ^ (-lA^Ag°xp . ..xA~i', 

TrGTlr- 


where 11^ is the set of permutations of length r and <j(tt) is the sign of tt. We deduce 
that d^ar{0) = 0 unless k = tt G 11^.. In the latter case, 

r —1 

d^aAO) = n V! = (-l)"^"^sf(r - 1), 

j=o 

where sf(m) = flidi *• i® ^ super-factorial. In terms of the Barnes G function, see 
(DLMF 2014, §5.17), we have sf(m) = G(rn 2). 

Consequently, by Leibniz’s formula, 


d^Aio) = E n 

r—1 


kl,i + ^2,; 
klA 


dtar{0)d^jar{0) 

E n { - 1 ), 

- h i = n V / 


T^l-\-TZ2 — k i = 0 


with 7ri,7r2 G 11^ both permutations of length r, i.e. the only terms surviving in 
Leibniz’s formula are those for which fci and ^2 are permutations. 

Let 

FAx) = ar{xy an-r{xY. (4.22) 
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Using the multi-index k = [fep,..., along with the definitions = [/cg,..., 

and = [fcr,..., fcn-i], we have 

dlF^[v) = 9^[a2(-l)a^,(+l)]=af’[a2(o)]af’[aL,(0)]. 

This derivative is nonzero only for -|- and where 

TTp' € and TTp' e n„_^. 

Combined with the above, we arrive at the expression 
d^Fa{v) = sf^(r — l)sf^(n — r — 1) 

,r'"47rP'=fe["l 

This expression is only semi-explicit and it is fairly difScult to proceed analytically 
with conditions of the form k = tti -|- 712 . However, the expression is valid for any 
fe, and we are only interested in the derivative that corresponds to the leading order 
term in expansion (14.171) . It turns out that, once all contributions are summed, such 
conditions drop out. 


(4.23) 



4.2 The leading order term 

4.2.1 Vertices with minimal weight 

We recall the definition of the function F in (14.151) . and we aim for the leading order 
term in (14.171) . which corresponds to the smallest m = \k\ such that d^F{v) does not 
vanish. 

It is clear from the preceding analysis that the number of (—l)s and (-l-l)s in v 
plays a crucial role. In particular, we call the difference (in absolute value) between 
the number of (-l-l)s and the number of (—l)s in a vertex the weight of that vertex. 
So far, we have been considering vertices with r (—l)s and n — r (l)s, so the weight 
is \n — 2r\. 

Note that for a derivative of order k = -|- -|- , where € H^ and 

[2l 

tt) G n„_r, we have 

|fc| = (r- — l)r -I- (n — r- — l)(n — r), (4.24) 

since for tt G H^ we have |7r| = iniiilllll ^ 

It is straightforward to verify that |fc| is minimal for vertices with minimal weight. 
This leads to r = n/2 for even n, and r = (n -I- l)/2 or r = (n — l)/2 for odd n. 
Consequently, we have to distinguish between these two cases. 
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4.2.2 The even case n = 2N 


Let us examine the factors in expansion (14.171) . There are 


vertices with r (—l)s, 


permutations of (I4.19|) . Let us call this set Vn,r- In this case, because r = n/2, we have 
r = n — r = N. For each vertex with minimal weight, i.e. for each v G V 2 n,n-, we have 
(—!)'’('") = (—1)'’ = (—1)^ and e'"^” ^ = 1. Furthermore, Pn,r{vY = 4’'("“’') = 4^ 
and / was assumed to be a symmetric function in its variables, hence / is constant 
on V 2 n,n- 

We have to sum over all derivatives of total order |fc| = 2N{N — 1). Since |fc| is 
minimal, we have from Leibniz’s formula that 

d^Fiv)=4^^fiv)d^F^iv), 

where 4 ^^ is the contribution of n- Each possible k is reached by a combination 
of permutations of length N. From ( 14.2311 . we find that 

= sf^(iV-l)sf^(Af-1) 


X E E n 




N-1 


X E E n 


J2]n 


2=0 

N-1 


2=0 


Identifying a sum with a determinant and permuting rows, 




7r2ens 2=0 


+ 7>'2,z 
'^lA 


= det(T[,'‘y^-)ij=o....,s-i 

= (-l)'^(-^)det(y:J.)..,=o......-i, 


where ^ j , L j = 0,..., s - 1. 

It is easy to see that det = 1. Indeed, it follows from the definition of that 
4W _ 4W — /iW 7 -i — 1 q - 1 

so subtracting the (j — l)st from the jth column we have 


det T[®1 = det 


1 


O' 


1 


Consequently, detAl®! = detT^®"^!, and by induction we have detA^®! = 1 for s > 1. 
Alternatively, this result follows from identifying Al®! as a classical Pascal matrix, see 
for instance (Edelman & Strang 2004). 
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Therefore, 


N-l ( [1] [1] 


.[1] 


(4.25) 


TT^^^ellAr 


2=0 


and the same holds for and instead of 7r[^^ and . Consequently, 

Y^diF^{v) = m\i\N-\). 


Assembling everything in formula (13.111) . the term corresponding to m = |fc| = 
2N{N — 1) in the asymptotic expansion becomes 


(-1)^4^' 


9 Ar \ A-i^ 

\mrsiiN-ir = ^siiN-i)\ 


{2N)\ (-ia;)2^' V ^ 

and therefore we arrive at the following result: 

Proposition 2 Let n = 2N be even, let v S V 2 n,n be a vertex with weight 0 and let 
f{x) he a symmetric function of its n arguments. If f{v) ^ 0, then 




l2N[f] = ^f{v)siiN-ir + Oic 


-2N^-1 


)■ 


4.2.3 The odd case n = 2N + 1 

Similar considerations hold in the odd case. In this case, we may have r = N and 
n — r = N+ loT viceversa, but these two cases are symmetric. They correspond to 
the sets V 2 Ar+i,Ar and V 2 Ar+i,Ar+i respectively. 

For V £ V2Ar+i,Ar, we have (—I)®!’') = (—1)'^ and = e‘“. For v £ V2Ar+i,Ar, 

we have (-!)«(’') = (-l)^+i and = e-^'^. Furthermore, P 2 N+i,Niv) = 

and m = |fc| = 27V^ in both cases. 

Proposition 3 Let n = 2N + 1 &e odd, let Vi £ V 2 Ar+i,Ar and V 2 £ V 2 Ar+i,Ar+i be 
two vertices with weight 1 and let f{x) be a symmetric function. If f{vi) ^ 0 and 
fiv 2 ) ^ 0, then 

^ 2 iv+i[/] = - ^2iv(^+i)+i 4^^'^^^^sf(iV - l)^sf(iV)2[e-/(ni) - e--/(^ 2 )] 

^(^(^-2Ar(JV+l)-2)^ 

If in addition f is an even function we have 

l2N+i[f] = - l)hmrf{v,)smoj + 

(4.26) 

Using the symmetric integral in (j^with / = 1, we have established in Proposition[2] 
that the Hankel determinant /i 2 Ar-i is strictly positive, for sufficiently large w. Hence, 
the next even degree polynomial p 2 N exists, because of (11.31) . Proposition [3] shows 
that h 2 N vanishes approximately at multiples of tt. Hence, the odd degree polynomial 
P2N+1 does not exist for these values of w. We examine this phenomenon in more 
detail in the next section. 
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5 The kissing pattern 


The kissing pattern of polynomials of consecutive degrees, as illustrated in Figure ??, 
results from degeneracy. The trajectories of the roots of an even-degree polynomial, 
say P 2 N-, and the next odd-degree polynomial P 2 N +1 tend to each other. These trajec¬ 
tories coincide at critical values of w, the roots of h 2 N- In that case, the lower-degree 
polynomial p 2 N satisfies the orthogonality conditions of an orthogonal polynomial of 
degree 2N + 1. 

5.1 Existence of the polynomials 

We recall formula (ira, that yields the monic nth degree orthogonal polynomial (x) 
in terms of determinants: 


Pnix) = 


^n — 1 


■ det 


Po 

pi 

f^n—1 

pi 

P2 

l^n 

pn 

Mn+1 ' * ' 

M2n-1 

ly at the zeros 

of hji—\ , 


1 

X 


(5.27) 


the analysis of the Hankel determinant. 

Using a different normalization we can define a polynomial that always exists, 
regardless of the zeros of the Hankel determinant. 


Pn{x) = hn-lPn{x). (5.28) 

This polynomial always exists but if hn-i = 0 then it has degree less than n. From 
the theory of quasi-orthogonal polynomials, or formal orthogonal polynomials, it is 
known that the degree of Pn equals the dimension of the largest leading non-singular 
principal submatrix of the Hankel matrix Hn-i (Brezinski 1980). This is equivalent 
to saying that the degree of Pn is equal to the degree of the first existing polynomial 
Pk of lower degree k < n. 


5.2 Degeneracy of P2 n+i 

We study what happens as uj tends to a critical value tJ. We assume that oj is suffi¬ 
ciently large, such that the even-degree polynomials exist but the odd-degree polyno¬ 
mial may not. Thus, we have /i 2 Ar(w) = 0, for some N > 0, and from (14.261) we see 
that this is a simple root. 

Recall the three-term recurrence relation (12.51) satisfied by the OPs: 

Pn+lix) = {X - an)Pnix) -/3nPn-l{x), U € Z+, (5.29) 

where p_i = 0 and Pq = 1. Using expressions (1^ we can rewrite this relation in 
terms of Hankel determinants, 

(^)^n—1 — i^hnhji — iX F i Pnix') ^nP^-^^x')- (5.30) 
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For critical values of w, the expression simplifies. Indeed, if = 0 and /i„_i ^ 0, 
the relation becomes 

h' 

Pn+i{x) =i--^Pn{x), (5.31) 

n.n-1 

i.e. Pn+i is a scalar multiple ofp„. This means that, at zeros of h 2 N, the polynomial 
P 2 N +1 reduces to a constant multiple of the polynomial p 2 N of lower degree. Hence, 
their zeros coincide and the trajectories of both polynomials ‘kiss’. The monic poly¬ 
nomial P 2 N +1 blows up, as its (27V -|- l)th zero - the one on the imaginary axis - 
necessarily becomes infinite. 

At this critical value w, p 2 N satisfies more orthogonality conditions than usual. 
The condition 

[ P27v(a;)p2Ar(a:)e'“"' dx = = 0 

J-l 'In-l 


implies that 


P2N{x)x^e^'^^ dx = 0 , 


'-1 


fc = 0, l,...,2iV. 


Thus, the polynomial xp 2 N still satisfies 2N orthogonality conditions. In the theory of 
quasi-orthogonal polynomials, it is customary to replace in this case p 2 N+i by xp 2 N, 
in order to obtain a complete basis for the space of polynomials. 

Such completeness is not necessary for the sake of quadrature rules and we forego 
a more complete description. Instead, we focus on an aspect of the kissing pattern 
that is unique to the polynomials at hand: the kisses in Figure ?? seemingly occur 
closer and closer to the endpoints ±1 as uj increases. This behaviour is studied next. 


6 Behaviour near the endpoints ±1 


The leading order term in the asymptotic expansion of the Hankel determinants /i„_i 
is very revealing, insofar as the existence of the polynomials is concerned. Yet, the 
expansion carries more information. It is conjectured and motivated in (Asheim et 
al. 2014) that the polynomials are approximately a multiple of Laguerre polynomials 
near the endpoints ±1 - recall expression (12.81) - and this is the property we investigate 
next using Hankel determinants. 

In other words, the observation is that for large ui the zeros of behave like 


±1 - 


iw 


where c is a zero of the Laguerre polynomial L|^„/ 2 J ■ In this section we intend to show 
that the leading order term in the expansion of (l — j^) indeed vanishes if and only 
if LLn/ 2 j(c) = 0. 

Using Heine’s formula (j4.18|) and the symmetric integral (14.1711 . we have 


-i n—i 

Pn = —Jn TT (l - - Xm) 

V iLuJ n\ \ luj J 


m=0 
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Thus, we invoke the theory of ^ again and proceed by expanding /n[/]i where 


f{xo, ■ ■ 


^n—1) — 


n—i 

)=n( 


m—O 


1 - - 

lUJ 


is again a symmetric function, but now depending on uj. 


6.1 The case n = 2N 

The analysis in section ij4] needs minor modifications in view of the fact that the 
function / itself now depends on w. The derivatives oi F = f have an expansion in 
inverse powers of u! and we have to take this into account. 

Recall that contributions to the asymptotic expansion (I3.11|) can be thought of as 
originating from vertices v S Vn,r with r (—l)s and n — r (+l)s. Our first question is: 
which vertices contribute to the leading order term in the expansion of /n[/]? Before, 
it was V G V 2 jv,jv- Here, in spite of the fact that / depends on oj, little changes. The 
leading order term still originates in vertices with minimal weight. 


6.1.1 The leading order term 

Consider a general vertex v G V 2 N,N+t, with 0 < t < N. Without loss of generality, 
and using the same multi-index notation as in the previous section, we can take v = 
( —1^^], Since / is linear in all its components, it is clear that d^f{v) = 0 unless 

ko,ki,... ,kn-i G {0,1}. Thus, suppose that G S such that 

kf\kf^ G {0,1}, = Ki and = K 2 - Then 




(-1) 


K1+K2 




V iuj) 

^ iujJ 


^_^'jKl+K22N + t — Kl — t — 


( iw )^-‘-«2 




(6.32) 


Note the absence of symmetry here: the roles of Ki and K 2 are not interchangeable 
because / focuses on the right endpoint a; = -1-1. 

There are three contributions to the leading order exponent of uj~^ in (13.111) : 


1. The dimension contributes n = 2N. 


2. The least order non-vanishing derivative of Fa, recall (14.221) . consists of permu¬ 
tations of length N + t and N — t respectively. From (14.241) . this contributes 
|A:| = 2(N — l)N + 2t^ to our exponent. 

3. A derivative of degree ki + K 2 contributes to the term 

in (j3.11|) and, from (I6.32() . additionally contributes - altogether 

We can choose ki = 0 (i.e., all derivatives can be only with 
respect to the trailing N — t components) and the contribution is N — t. 

The total exponent is 2N‘^ -|- -|- 2t^ — t and this, clearly, is minimised for t = 0. 

Though derivatives of / may contribute positive powers of uj, vertices with larger 
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weight (i.e., larger t) contribute smaller powers of uj, and the latter effect is stronger. 
The resulting leading order behaviour, , is a factor smaller than that in 

Proposition [2] simply because that is the size of / at an endpoint with weight 0. 


6.1.2 Analysis of the derivatives of F = fg^ 

Let us examine the derivatives of F = fg^ further. Corresponding to a vertex v G 
V 2 N.Nj the leading order term must differentiate g^ with permutations of length N. 
In addition, we may have derivatives of / with respect to the second set of variables. 
Denote by G {0,1}^ any vector such that |Ps| = s, for s G {0,..., A}, then we 
compute the derivatives from the above formula: 


I' — l ( — 

^ - (1 + Oiu-^)) 


(iw) 


(6.33) 


which corresponds to the choice t = 0, ki = 0 and K 2 = s in (16.321) . Derivatives 
of the former we have already analysed in From (I4.23[) . and remembering that 

so we find that 



22^%f4(iV - 

1) 




N-1 

X 

M 

_^^(T(7ri) + CT(Tr2) 

n 






N-1 

X 

M 

_]^^0'(7r3)+0-(Tr4) 

n 

7r3+7r4=fct^^ 




'^ 1,1 + '^ 2,1 
'^ 1,1 


■^ 3,1 + '^ 4,1 


Note that the first sum in the previous formula is equal to A!, because of (14.25^ . 

The increase of the order of the derivative of F in expansion (I3.11|) comes at a cost 
of a factor (—iw)”®. On the other hand, a higher order derivative of / yields a factor 
of (ia;)“^+®, from (j6.33|) . The product of these factors is (—l)®(ia;)“^ and has the 
same asymptotic size in w for all s. Hence, we need to consider all s, 0 < s < A. 

Since derivatives of g^ vanish unless the order of the derivative is a combination of 
permutations of length A, we need to consider all such combinations of permutations 
and all values of s. This leads to a sum of terms of the form 




Ps,i = l 


for s G {0,1,..., A} and with and sums of two permutations in Hv- 
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The sth term in this large sum is 

E E n +1)3?"’ 

|p^|=S 7ri,7r2,7r37r4£nN' Pa,i = l 

= 2^^^+^si\N-i){-iy(ic,r-^{-cf-^ E n + 

|p^ |=s 7^1.7^2.7^37^4G^^;• Ps,i—1 

N —1 y \ TV — 1 y 

_l',tT(7ri)+t7(7r2) JJ ( 7ri,i+7r2,i \ (_^-^<T{7T3)-\-<Ti7T4) 77 ( ^3,i+7r4,i 


2 = 0 


s-N/ \N- 


i-iy 


n 


i=0 


T^3,i 


'sr(iV - l)(-l)"(iw)"-^V-c) iV! 

7V-1 


X(-1) 

^ 22N^+N^f4 

X ^ ^ (-iy(^3)+cr{7Vi) JJ- (^3,i + TT4,i+Ps,iy- 

Ip^I^s 7r3,7r4£nK *—0 

In the last computation, we have used (14.251) and the fact that 


N-l 

cri'!r3)+<yi'rri) I I ( '^3,i + '^4,i 


n 

2=0 


E n (’^3-*+ ^4,i+!)(-!) 

•7r3,Tr4Gn;v Pe,i = l 

^ ^ (_l)<T(773) + <T(7r4) JJ- (^3,^ + 7r4,i +Ps,i)! 

7r3,7r4€niv 


7^3,2 


2=0 


^3,2l7r4,^! 


In this last sum, every consists of s ones and N — s zeros, hence there are (^) 
such vectors. Note that each gives exactly the same result, because everything else 
in the relevant expression is constructed from two permutations. Therefore, we might 
just consider 


N—s times s times 

p: = ( 


/N 


) times. Therefore 


N-l 


^ ^ (_l)^(7r3)+<7(7r4) JJ (^3,2 +7r4,2 +Ps,0! 

Ip^I^s 7r3,7r4eniv i—0 

^ '' 7r3,7r4enN ’ ’ 


(-l)'"^’^Met£’[^’"l(7r), 

Trelliv 


where 




(7r» + j)! 
Trjj! 

(TTj + j + 1)! 

7rdj! 


i = 0,..., — s — I, 




j = 0,...,N-l. (6.34) 
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Proposition 4 Let tt = with 

= (tto, ... jTTat-s-i), Trt^J = (ttat-s, ... ,7rAr-i). (6.35) 

If ttI^I G Htv-s and tt^^I G — s + Itg, then 

N\‘^ 

= (6.36) 


otherwise det f ( tt ) = 0. 

Proof Suppose first that the hypothesis on tt^^I and holds. We deduce, 
rearranging rows separately in the first N — s and the last s rows of f [^’®l(7r), that 

(_ 1 )'^(^) det £■[^’'’1 ( tt ) = det (0, 1 ,..., - 1 ). 


Since 


c[N ,s] 


(0,l,...,iV-l) 


|C:o 


i = 0,..., N — s — 1, 
i = N — s,..., N — 1, 


and j = 0,..., N — 1, once we extract a factor of i + 1 from rows i = N — s,..., N — 1, 
the outcome is 


- 1 ) 


Nl 

{N-s)l 




where 




( 0 , 1 ,.. 





V *+i y’ 


i = 0,..., A^ — s — 1, 

j = 0,...,iV-l. 

i = N — s,..., N — 1, 


The matrix C is somewhat easier to manipulate. First, note that equals the 

Pascal matrix used in the analysis in 114.2.21 hence detC^^’^l = 1. 

In case s = 1, easy calculation with binomial numbers confirms that 


and 


Therefore 


j j — 1 — ^i— 1 j ’ 


j = 0,...,iV-2 


^IN,1] _^IN] 






detCl^’^l =det 


1 0^ 
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where both matrices in the lower right block differ in only one column. This leads to 

detCl^’^1 = detC[^-i’^] +detC[^-i’°l = detCl^-^l + 1 


and we deduce that detC^^’^l = N. 

Let us generalize the above to larger s. Note that 




'^[AT—l,s] ^[Ar,s —1] 

^[Ar-l,s] [N-l,s-l] 

^[Ar-l,s] _ ^[Ar,s-1] 


j = 0,... ,iV - s - 1, 
j = N -s, 

j = N-s + l,...,N-l, 


An argument identical to the one we have used before shows that 

detC'^’"' = detC'^-^’"' + detCl^"^’""^!, s = 0,1,..., A^. 


In case N = s, the same reasoning leads to detC^^’^l = detC^^’®! = 1. Induction then 
shows that detCl^”'’! = (^), and (I6.36P must be true. 

Suppose now that tt is not isomorphic to IItv-s © (N — s + IIs), which makes sense 
only when s € {I,...,iV — I}. Then there exists (at least) one integer r > N — s such 
that TTr G and = iV — s — i for some i >1. Among those, we choose r so that 
TTr is minimum, and we take t to be the index such that iTt = + 1. Then either 

TTt & or TTt G In the first case, using (16.341) . the r-th row of £'[^’®l(7r) is 


(A^ — s — z + j + I)! 
(N — s — i)!j! 


and the f-th row of is 


(A^ — s — z + j + 1)! 
(Af-s-f + l)!j! ’ 


j = 0,...,iV-l, 

j = 0,...,iV-l. 


Since the t-th row is a scalar multiple of the r-th row, the determinant vanishes and 
the proposition is true. If zr^ G then we repeat the previous reasoning with the 
index t and the index u such that 7r„ = tt^ + 1. We continue this process and at some 
point we must find two indices with the property above, since it cannot happen that 
all permutations in tt^^I take smaller values than those in □ 


6.1.3 Asymptotic behaviour of p 2 N 

Let us assemble everything. Since there are precisely s!(JV — s)! permutations iso¬ 
morphic to IIjv-s ® (A^ — s + IIs) (hence, which produce nonzero determinants in the 
proposition), we obtain 


(_l)'^(^3)+<^(7r4) JJ" ('^3,i+'^4,i +Ps,i)'- 


|Pg|=S 7r3,7r4€n2V 


2=0 


’^3,2*^4,2- 


s!(A-s)!2 
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for s = 0,..., iV. 

The contribution of v includes the derivatives derived above, as well as an addi¬ 
tional factor (—iw)“® for each s arising from the 0 ( 0 ;“™“") term in (13.lip , and the 
factor (—= (—1)^. This totals, after some manipulations. 


N 


\N-s 




^(_l)^(_i^)-22^ +^sf(iV - l)4(-l)^(ia.)-^(-c)^ _ 

s—0 '' 

= (-icc)-^22^'+^sf(iV - 

e=0 ® / 

,sf(iV)4 


= (-iw)-^2 


-Nr)2N^+N^ 


m 


-LAr(c), 


where the latter simplification follows from the known explicit expression 



(6.37) 


See for instance (Chihara 1978, Chapter V) or (Szegd 1939, Chapter 5). 

There are vertices in V 2 N,n, hence the leading term in the expansion of 

P 2 Af(l - c/(ia;)) is 

22iv^+iVsf(jv)4 ^ , , , .,^ 2^^^+^sf(7V)% 

(-ia;)2^"(-ia;)^iV!3^^^^ ^ N\^ 

Finally, we can divide by the leading term of h 2 N-i, recall Proposition [21 

4^'sf(Af- 1)4 

to obtain the leading term of the monic polynomial: 

P27v(l-—)=- jv—^LAf(c) (1+ 0 ( 0 ;“^)) , UJ^OO. 

\ lUJ/ ^ ^ 

This is precisely the result we wanted to show. 


6.2 The case n = 2N + 1 

Now, we consider a general vertex v S V 2 N+i,N+t- As before, d^f{v) = 0 unless the 
multi-indices are such that £ {0, 1 }'^“'-*^ fcPl g |q^ with = ki and 

= K 2 . In this case, we have 


, , , , /_1 '\Ki-I-K 2 oAT-l-t-Ki/_ 'lAr-|-l-t-K2 

e'"m = ‘ -(>+^ 


The contributions to the leading term are as follows: 


(6.38) 


1. The dimension contributes 2N + 1. 
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2. The least order non-vanishing derivative of fa, that consists of permutations of 
length N + t and N + 1 — t. This amounts to |fe| = 2N'^ + 2t{t — 1). 

3. By similar reason as for n = 2fV, from (I6.38L the leading power is ki -|- K 2 (the 
degree of the derivative) plus N +l — t — K 2 (the above contribution) - altogether 

-|- 1 — t -|- Ki- Since we wish to minimise this, we need to take ki = 0 and the 
contribution is N — t + 1. 

The total exponent is therefore 2N‘^ + 3N -1-2-1- t{2t — 3), which is minimised for 
t € {—N, ..., A^-|-1} in the case t = 1. The exponent is then equal to {N + 1)(2A^-|-1). 
Therefore, we need consider just vertices in V 2 N+i,n+i- 

As before, derivatives of / in the term F = fg^ may have order s > 0. It is 
important to observe that the range of s is unchanged, i.e. 0 < s < A^, since it 
is constrained by the number of (+l)s in the vertex v S V 2 Af+i.Ar-i-i which is N as 
before. Going through similar computations, the identity (I6.37|) quickly resurfaces in 
the leading order term. 

This leads to the following result. 

Theorem 5 For uj ^ 1 the zeros of P 27 v(') + ^(^~^); where 

..., > 0 are zeros of the Nth Laguerre polynomial. The same is true for 

P 2 n+i> except for a single zero on the pure imaginary axis. 

Proof We have shown the asymptotic agreement with the roots of the Laguerre 
polynomial. The single zero on the pure imaginary axis follows at once from the 
symmetry of with respect to the pure imaginary axis. □ 


7 Roots of hn in the complex plane 

In this section we investigate in more detail the zeros of the Hankel determinants in 
the complex plane. First of all, we recall that h„(aj) is real for a; € R, so all complex 
zeros must come in conjugate pairs. Also, since /i„(—w) = h„(a;), we can restrict 
ourselves to the first quadrant of the complex plane. 

Figures 17.31 and 17.41 show these zeros for different values of n. They follow very 
regular and symmetric patterns reminiscent of onion peels, that we intend to explain 
in this section, at least for large values of uj. These patterns result from a delicate 
balancing act, in which algebraic powers of become comparable in size to decaying 
complex exponentials of the form e“^'^^, c > 0, in the upper half of the complex plane. 
We revisit the asymptotic analysis of (21 this time taking complex exponential factors 
into into account. 

As before, we need to distinguish two cases, corresponding to even and odd values 
of n. 

7.1 The odd case: roots of h 2 N-i 

We commence from h 2 N-i = l 2 Ar[l]. The contribution to h 2 N-i in the asymptotic 
expansion (13.111) corresponds to the layers V 2 N,N±t, t = 0,... ,N, with the previous 
notation. We quantify the contribution of each layer with two numbers: 
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Figure 7.3: The zeros of hi, h^, and in the complex plane. 





Figure 7.4: The zeros of h 2 (w), h 4 (u;) and h(,{uj) in the complex plane. 
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a) a complex exponential factor: each vertex v G V 2 N,N±t contributes = 

b) the leading power of ui ‘originating’ in vertices in the layer. As computed before, 
recall (I4.24|) . that contributes a total power 2N‘^ — 2N + Therefore, in the 
end we obtain ^- 2 N- 2 iN^-N+t^) ^ ^- 2 N^- 2 t\ 

In other words, each layer V 2 N,N±t contributes 

P2nM^) ■■= ^2NW + ^(^■')] • ( 7 - 39 ) 

We now need to choose a plus or a minus sign: e“^ is large for Imw < 0 and e““^ for 
Imw > 0 , so we choose the latter, since our goal is to balance it with a positive power 
of uj. We thus assume that e”’*^ = 0{ujP) for some p > 0. The immediate implication 
is that 

P2nA^) = 

The idea is now to choose ti < t 2 in {0,1,..., iV} such that P 2 NM and P 2 N,t 2 are of the 
same order of magnitude and the remaining P 2 N,t are of smaller order of magnitude. 
The first requirement means that 

-2N‘^ - 2tl + 2tip = -2N‘^ - 2t\ + 2 t 2 P, 

therefore p = + ^ 2 , and then 

The second requirement is that 

-2N'^ - 2A + 2t{ti + t2) < -2N‘^ + 2tit2, 

therefore {t — ti){t — t 2 ) >0 for all t ^ ti,t 2 - 

Let ti = k, t 2 = k + 1, where /c = 0,1,..., iV — 1. Then p = 2fc + 1, 

P2N,k,P2N,k+l = 

P 2 ^.t = o(a;-2^'-2*'+4‘'=+2*), + 

It is trivial to verify that indeed {t — k){t — k — 1) > 0 for all t 7 ^ fc,/c + I, so all 
is well. Moreover, these are all possible such choices. For suppose that there exist 
ti,t 2 such that ti < t 2 — 1. Then (t — ti){t — ^ 2 ) < 0 for t = ti + 1 and we reach a 
contradiction. Thus, we have exactly the right number of N different choices, which 
correspond exactly to a separate ‘onion peel’ in the lower-right quadrant in the Figure 

[7:^ 

Thus, we choose k G {0,1,..., A^ — 1} and extract just the fcth and {k + I)st terms, 

g— 2 icj/c g— 2 icj(fc-l-l) 

h2N-lAA = C27V,fe^2Ar2+2fc2 + ,k+l 

^-2\u,k / g- 2 ia; \ 

“ ^^ 2 Af 2 + 2 fc 2 y^2N,k + ^2N,k+l ^2{2k+l) J ' 
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Setting h 2 N-i,k{^) to zero, we obtain an asymptotic expression for the fcth ‘onion 
peel’. This leads to the equation 


^2{2k+l) ^2iuj ^ C2N,k+l 
C2N,k 

We compute the coefficients CN,k later on, in order to arrive at explicit expressions 
and rule out the case C 2 N,k = 0. 

Taking roots, 

/ \ l/(4fc+2) 

^gW(2fe+l) _ I- 2_N^\ ^nie/{2k+l)^ e = 0,l,...,4k + l. 

\ C2N,k ) 


Therefore 


^ia;/(2fc+l) 


C 2 Af,fe+l 

2/c + 1 2/c + 1 C2N,k 

and we deduce that 


l/(4fc+2) 


^7rii/i2k+l) 


OJ 


-(2fc+l)iW ^ 

y 2A: + 1 \ C2N,k J 


X l/(4fe+2) 


7rif/(2fe+l) 




= 0,l,...,4fc + l 


£ = 0,l,...,4fc + l. 


(7.40) 

in terms of the Lambert W function, see for instance (Olver, Lozier, Boisvert & Clark 
2010, §4.13). 


7.2 The even case: roots of h 2 N 

We revisit the work of the last subsection, and now we have the layers V 2 N+i,N-t and 
V 2 N+i,N+i+t, t = 0,..., N. The exponential factor is in the first case, and 

g-( 2 t+i)n.j second one. 

For reasons of symmetry, it is enough to look at one of these cases. Since we 
are interested in the upper-right quadrant, we choose the second case. Computing as 
before, the leading power of w is 21V^ -|- 2t^ + 2t. Adding the dimension n = 2N + 1, 
we deduce that the contribution of V 2 N+i.N+i+t is 

g—(2i+l)icA; 

P2N+iA‘^) := ^2N^+2N+l+2f^+2t [c 2 N+i,t + 0{ui~^)] , r = 0,..., - 1. (7.41) 

Assuming again that e““^ = 0{uj^) for some p > 0, we have 

P2N+l,N+l+t = 0(^-2V=-2iV-l-2t=-2t+2pt+p)^ 

We again choose ti < t 2 according to the two rules. This leads to p = -I- ^2 + 1, and 

p p _ coo 1 — ‘2N^ — 2V-t-t 1-t-t2-t-2ii^2 “v 

P2N+i.N+i+t = 
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In particular, the second requirement reduces again to (t — ti)(t — ^ 2 ) > 0 for all 
t 7 ^ ^ 1 :^ 2 - It follows at once that ti = fc, ^2 = ^ + 1 for some fc € {0,1,..., iV — 1}, 
otherwise the inequality fails for + 1. 

So, all that remains is to check if the above choice of consecutive ti,t 2 works and 
indeed, trivially, it does: either t < ti ,^2 or t > ti ,^2 and in each case the inequality 
works. In other words, p = 2k + 2 and we investigate the zeros of 

g—(2fe+l)iLi; g—(2fc+3)itiJ 

h27V+l,fc(w) — C2iV+l,fc^2Ar2+2Ar+l+2fe2+2fc + ^2N‘^+2N+l+2k^+6k+i 

g-(2fe+l)iiii / \ 

~ ^2N^+2N+l+2k^+2k ( 02Ar+l,fc + C2Ar+i,fc+l ^4j,_|_4 j ■ 


Therefore 


^4A:+4g2ia; ^ _ 02A:^+l,fc+l 

C2Ar+l,fe 


iuj / {2k-\-2) 


Lue 


and we deduce that 


uj Rs -2i(fc + 1)W 


C2Af+l.fc+l 
C2Ar+l,fe / 


C2Ar+l,fe+l 


X l/(4fc+4) 


^nil/(2k+2)^ ^ = 0,...,4fc + 3 


l/(4fc+4) 


^TTie/{2k+2) 


)■ 


. , V. ., ^ — 0,..., 4A: + 3, 

2(fc+l) \ C2Ar+i,fe / / 

(7.42) 

again in terms of solutions of the Lambert W function. All this explains, at least in 
an asymptotic sense, the ‘onion peel’ structure of zeros of /i„. 


7.3 The coefficients ^ 

The coefficients Cn,k are defined in (I7.39P and (17.411) as the contribution of the layers 
of vertices V 2 Ar,Ar+fc or V 2 Ar+i, 7 V+i+fe- 

Elaborate calculations, similar to the ones in (JU lead to the following results. For 
odd n = 2N — 1, it is the case that 

C 2 N,k = 4^'-'='[sf(iV -k- l)si{N + k- l)]^ 


and consequently, see (I7.40L 


/ „ \ l/(4fc+2) „ 7 ri/( 4 fc+ 2 ) 

/ C2N,k+l \ e ^ ' 

■ (iV + k)'. ' 

V C2N,k J 2 

_iN-k- 1)!_ 


fc = 0,...,A^-l. 


In the even case, we have 

C2Af+i,fc = i(— ^■^N+k^{N-k)iN+k+l) fc)sf (iV -k-l)f 


consequently, see (j7.42|) . 


/ C 2 Ar+l,fe+l^ 

l/(4fc+4) g7ri/(4fc+4) 

■(Af + fc + 1)!' 

\ C2N+l,k / 

2 

{N-k-iy._ 
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Figure 8.5: On the left, plot in log-scale of |hi(a;)| (solid line), 1 ^ 3 ( 0 ;)! (dotted), |h 5 (a;)| 
(dashed) and 1 ^ 7 ( 0 ;)] (dashed-dotted). On the right, plot in log-scale of |h 2 (a;)| (solid 
line), 1 ^ 4 ( 0 ;)! (dotted), |h 6 (a;)| (dashed) and |h 8 (a;)| (dashed-dotted). 


8 Roots of hn on the real line 

In this section we study the zeros of the Hankel determinants for w > 0. From 
the information contained in Figures 17.31 and 17.41 one can expect different results 
depending on the parity of n. Indeed, if we compare different values of n, we obtain a 
picture like the one shown in Figure 15751 where we show the value of h„ in logarithmic 
scale as a function of uj. 

Based on these results, and similar ones that can be obtained by direct computation 
in Maple, we formulate two results about the properties of the Hankel determinants: 

1. Two consecutive Hankel determinants h„ and hn+i cannot have any common 
positive real zeros for n > 0 and w > 0 . 

2. For n > 0 and w > 0, the determinants hn and h „+2 do not have any common 
positive real zeros. 

We also present two conjectures: 

1. For n > 0 and w > 0, the determinant h 2 n+i does not have any positive real 
zeros. 

2. For n > 0 and w > 0, the determinants h 2 n and h 2 n +2 have interlacing zeros. 

Note that the first conjecture implies that P 2 n exists for all w > 0. 

An auxiliary result is the following important identity, which is known (cf. for 
instance (Bleher & Its 2005, Section 2)) but is presented for the sake of completeness 
in the next lemma. 
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Lemma 6 It is true that 


(^n) — ^ ^ 1; 


(8.43) 


where ' indicates derivative with respect to ui. 

Proof This identity is closely related to the Toda differential equation, and it is 
of great importance in integrable systems and random matrix theory, see for instance 
(Bleher & Its 2005, Section 2), (Bleher 2011, Proposition 18.1) or (Bleher & Liechty 
2013, Theorem 1.4.2). For completeness, we present here a brief proof in the current 
setting, which uses the connection between the Hankel determinant and the partition 
function of a certain random matrix ensemble. 

We revisit the Heine formula for h„_i, see (13.9L and we note that using elementary 
operations, we can replace the row (xq, xl,..., for i = 1, 2,..., n — 1 in the 

Vandermonde determinants by (pi(xo),Pi(xi),... ,pi(x„-i)), provided that all Pi(x) 
exist. This leads to an alternative formula for the Hankel determinant: 


n —1 

hn—l — , 

3=0 


Ki = 



iujx 


dx. 


Next, we observe that 


xpj{x) 


2^iujx 


dx = 


Xjrvj, 


using the recurrence relation (12.5|) and orthogonality. As a consequence. 


(8.44) 


(8.45) 


n — 1 

^n-l = Y. 



n-l / 

= E n' 

e=o \j^£ 


n — 1 

iaiKi = \hn-i ^£- 


(8.46) 


We differentiate again, bearing in mind that = i(/3f+i — (it), see (12.71) : 


n—1 n—1 

= i/i(,_i Y o^£-hn-i ^(/3£+i-/3£) 
(.=0 £=0 


{K-lf ^ ^ i.K-1? Khn-2 

— - tln-lPn — —r -T- 

lln-l On-l lln-1 


where we have used the second equation in (12.61) and telescoped the second sum, taking 
Po = 0. Multiplying throughout by hn-i and shifting the index n — 1 i—n, we obtain 
the result. □ 


We note that (18.431) can be rewritten as 

d^ 

l) “ 

in terms of the recurrence coefficient /3„. 

We first show that no two consecutive Hankel determinants can vanish simultane¬ 
ously. 
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Proposition 7 There is no n>l and w* > 0 such that 

Proof Assume that = hn{u)*) = 0 for some value w = w*. Then by (18.431) 

we have = 0 and so has a double root. Since both terms on the left hand side 
of (18.431) have a double root, so must the right hand side and this implies that either 
hn+i = 0 as well, or that /i„_i has a double root. 

In the latter case, two consecutive Hankel determinants have a double root. In the 
former case this happens too. Indeed, in this case we have hn = h'^ = /i„_i = = 

0. We can reformulate ()8.43|) as 

h^j^-^hji+l ihjij^i) — hnhn^2: nPO. 

It follows that h'^_^_^ = 0, i.e. both hn and hn+i have a double root. 

It remains to rule out two consecutive double roots. Let’s assume they are hn and 
hn+i- In that case, the right hand side of (15.301) vanishes at w = w* but the left hand 
side does not, since does not vanish identically, unless also h„_i = 0. Subsequently, 
we can deduce from another reformulation of (|5.30l) that h'^_i = 0 too. Continuing 
this reasoning leads to a chain of double roots and all Hankel determinants vanishing 
down to n = 0, which leads to contradiction. □ 

It follows immediately that one can not have hn = h'n = 0 either, since by (18.431) 
at least one of hn-i or hn+i has to vanish too. 

We can also exclude hn = hn +2 = 0, which is our second result: 

Proposition 8 There is no n > 0 and to* > 0 such that 


hn{u}*) = hn+2{uJ*) = 0 . 

Proof The result is true by direct computation for n = 0 and n = 1. Let us assume 
it is true up to n — 1, and assume that = 0 for some w* >0. We intend to 

show that /i„+ 2 (w*) 0. 

We know that hn-i ^ 0 by Lemma[7]and that hn -2 0 by our inductive assump¬ 
tion. It follows from (12.61) that an-i is analytic at w*. It also follows from (12.61) that 
j3n{oj*) = 0. Since h'n^O and hn -2 0, this root of /3„ is simple. 

We reformulate the differential-difference equations (IQ) as 


/^n+l 


Q^n+1 


—ia^ + Pn 

.P'r. 


n+1 


O-ri 


Plugging in a Taylor series of On-i and /3„ around to = uj* and using the above 
recusions shows, after straightforward computation, that an and an+i have a simple 
pole, Pn+i has a double pole, Pn+2 has a simple root and an+2 is analytic at uj* . Using 
the expressions 


_ {xPn,Pn) 
{Pn,Pn) 


and 


{Pn,Pn) 
{Pn—1 ; Pn — l) 
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this implies that {pn+i^Pn+i) has a simple pole and {Pn+ 2 :Pn+ 2 ) 7^ 0- The latter in 
turn implies that /1„+2(W*) ^0. □ 


We observe that this type of ideas were used in a similar problem (but with a 
complex cubic potential on a union of infinite contours in C) in the thesis of N. Lejon 
(Lejon 2016); in that setting, the combination of analogous properties for the Hankel 
determinants and the string (or Freud) equations lead to the proof of existence of the 
orthogonal polynomials of even degree. In the present case, main issue is that the 
string equations have a more involved structure, because of the presence of boundary 
terms. This topic is under current investigation. 
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